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Abstract 

The optical conductivity for the surface excitations for a Topological Insulator as a function of the 
chemical potential and disorder is considered. Due to the time reversal symmetry the chiral metallic 
surface states are protected against disorder. This allow to use the averaged single particle Green's 
function to compute the optical conductivity. We compute the conductivity in the limit of a finite 
disorder. We find that the conductivity as a function of the chemical potential /i and frequency O 
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is given by the universal value a{Vl > 2/i) = For frequencies Q, < fi and elastic mean free path 
lei = VT which obey kpl > 1 we obtain the conductivity is given by a{Vt < 2\iJ,\) = {Qt^^^+i • 
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the limit of zero disorder we find a{fi ^ 0, il, kpl — )• oo) = 



1 



I. Introduction 

For time reversal invariant systems one finds that for Kramer's states tlie time reversal 
operator i3 obey = — 1 and thus the second Chern number is given by (— 1)*^ = —1 where 
1/ is an odd number. Topological Insulators (TI) [l-4| are characterized by chiralic gapless 
electronic spectrum. For the 3D TI material Bi2Se3 the surface state consist of a Weyl 
equation with single Dirac Cone which is below the chemical potential fi and the bulk gap 
0.3eV. Due to the topology the backscattering is suppressed Q and therefore localization 
might be prohibited. Tunneling scanning has confirmed the presence of backscattering and 
therefore the presence of a topological metallic state. The conductivity results are less 



conclusive, due the presence of the 3D bulk gap {g] or the insulating gap for thin layers TI. 
The charge current operator for Weyl equation is identified with the spin on the surface of 
the TI. Therefore for clean systems conservation of momentum will cause the low frequency 
conductivity to vanish, only the optical conductivity for frequency Q > 2fi {fi is the chemical 
potential) is finite. In the presence of impurities elastic scattering conserve energy but not 
the quasi momentum, giving rise to a finite Drude conductivity. Due to the spin orbit 
interaction anti-localization effects dominate the physics giving rise to a metallic conductivity 
jsl, liol 12, 18|. For the Weyl equation the velocity operator is given by: ^ = vpeija^ and 



= 1,2 and the charge current is given by the spin operator. We formulate the optical 
conductivity for a finite chemical potential. We define Green's functions for particles and 
anti-particles (holes). Our results dbTQ clS clS following: For a finite chemical potential we 
find a{il > 2/i) = ^[/^ + ^.fp.D.i—^ + /^)] which at low temperature takes the universal 
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value a{uj > fi) = The results obtained here in the absence of disorder are similar to 



the one obtained in graphene 16|]. For weak disorder and chemical potential /i > we 
find the that the conductivity is given by the Drude conductivity, a{il < 2\fi\) = (n-r)Hi " 
{lei is the elastic mean free path). The plan of this paper is as followings: In section // 
we introduce the model. In section III we construct the Green's function for a model 
with a finite chemical potential. In chapter IV we consider the effect of the white noise 
scattering potential on the single particles Green's functions. In chapter V we compute the 
conductivity in the single particle approximation. Chapter VI is devoted to conclusions. 

II- The Weyl model 
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The surface state Hamiltonian of a Topological Insulator (T/)of the 8128 es family ma- 
terials is given by a the Weyl model [? ] . The presence a random potential and an electro- 
magnetic field modify the Weyl model in the following way: 

S = j dt[j (fr-^\f,t){idt- eAo{f,t)-^{f,t) - H]; 

J n n nv 

(1) 

V is the Fermi velocity ,A{t) — ^^(Q)e~*^* is the external vector potential, Ao{f,t) is the 
external scalar potential and Usd^r) is the random potential controlled by the white noise 
correlation function << Usc{k)Usc{k') >>= DscS'^[k + k']. In the absence of the random 
potential and external fields the Hamiltonian (in the first quantization) in the momentum 
space is given by h[k] — —a'^ki + a^k2- The Hamiltonian h[k] is time reversal invariant 
hlkl'&hl—k] ( •& — ia'^K is the time reversal operator and K is the conjugation operator 
which obey ■j?^ = — 1 ) This Hamiltonian has two eigenvectors : \u^'^\k) > with the eigenvalue 
Ep = = hv\k\ for particles and {u'^'^k) > for anti-particle Ea — — —hv\k\. 

\u^-){k) >= ^\k > ®[-l,ie'^^''^''yY-,x{kx,ky) = tan-' (2) 

The spinor operator ^(r) is decompose into the eigen modes of the unperturbed Weyl 
Hamiltonian. 

^(f) = ^e*"lC(fc)M(+)(fc) + B\-k)u^-\-k)] =J2e''''"^{k) (3) 

k k 

H = hvkp stands for the chemical potential ,|0 >= > ®|/u;a >. > stands for 

particles ground state and obeys C{k)\fi;p >= for E(k) > /j, = fivkp > {kp is the Fermi 
momentum. For < we have in the ground state antiparticle for energies E(k) < \fj,\). 
\fj,; a > represent the antiparticles ground state and obeys B{k)\fj,; a >— 0. Using the eigen 
spinors the Hamiltonian is equivalent to two coupled bands. 

Ho = Y^[hv\k\C\k)C{k) + hv\k\B\-k)B{-k)] 



Hd^Y^Y. Usc{k - P) [VP'nk,p)C^{k)Ci^ + V^'^-k, -^B{-k)B\-p) 

k P 

-^C\k)B\-p) + l/«'f(-fc,p)S(-^)C(p)] 



(4) 



Ho is the unperturbed Hamiltonian and Hn is the effect of the random potential on the two 
bands. The spinors structure gives rise to vertex functions for the coupling to the random 
potential in momentum space The vertex functions are given in terms of the particles and 
anti-particles matrix elements where p stands for particles and a stands for antiparticles: 
are given by: 

VP'P{k,p) =< u^+\k)\u^+\p) >; V'^'i-k, -p) =< u^-\-k))\u^-\-p) >; 
VP'''{k,-p) =< u^+\k)\u'^-\-p) >;V'''P{-k,p) =< u^-\-k)\u'^+\p) > 

(5) 

In agreement with the time reversal invariance -i?^ = — 1 we find that the backscattering is 
prohibited: VP'P{k,p = —k) = e^^^^^)^^^^'^)) cos[|(x(^) ^ xi^k))]. From the representation 
x(^i)^2) = x(^) = tcLn^^{^) we obtain the relation x{~k) = vr + and find that 
VP'P{k,p= -k) = 0. 

III-The single particle Green's function for a finite chemical potential yu in the 
absence of disorder 

In the absence of disorder we can use the spinor operator \E'o-(^) = [^a=^{k),'^^=i{k)'^] 
with the components : 

^^{k) = C{k)u^+\k) + B\-k)u^-\-k) (6) 

Using the spin half spinor \a = ±|) we obtain the the projected spinors u^^\k) and u^~^(k) 

ul+\k) = {a\u^^\k)), ui-\k) = {a\u^+\k)) (7) 

We define the single particle Green's matrix igcr,cr'{k, k'; t,t') = i6^^,Gf^^„i{k] t — t' = r) [isl - 

tG^,^,ik,T) =< 0|T(vl/,(^,r)vl/t,(^,o))|0 >= iul^\k)r ■ uiVik) < t,;p\TiC{k,r)C\k,0))\f^;p > + 
{ui-\-k)y -ul^^-k) < fi;a\T{B\-k,0)B{-k,T))\i2;a>= 

(8) 

The Fourier transform in the frequency domain allows to write: G^{k, u = 00+ fi) for particles 
and G"'{—k, u = u+fi) for antiparticles ( p stand for particles and a for anti-particles Green's 
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function). The Green's functions for a positive chemical potential > are given by in 
terms of the Fermi Dirac function for temperatures T — >■ 0, fp.D.ii^) — -^^i^^ii', /5 — >■ oo 

ryprt , n^ ^ - fp.D.i^) , fF.D.{uj) ^„ 
G^ik.ui — cu + a; a > 0) = 1 , e — >■ 

u-{E{k) - ii)) + ie {E{k)-n)-ie 



ui + hj[k) + n — le 



(9) 



IV-The single peirticle Green's function in the presence of the white noise 
potential Usc{r) 

To second order in the scattering potential Usc{r) the self energy for the particles is given 
by T,{k,u;p) for the anti-particles the self energy T,{k,u;a) is given by a similar result: 

E(^,a;;p) = J ^ « Usc{k - k')U,,{k' - k) » VP^^{k,k')V^'^{k\k)G^\k' ,oj) 

+ / ^ « Usc{k-k')Usc{k'-k) » VP'''{k,-k')V''^P{-k',k)G^''\-k\u) 

u; p) = T,R{k, u; p) + iLi{k, u; p) 
^n{Ku;p^u + t.,A;p)^r,j^ -^-^^ + . ^ ^ 

T^i{k,u;p) = -7rr|(a; + n)[fF.D.{-^) - /f.d.(^)] - 7rr^((j + fi)fF.D.{-uJ - fi) 

n - ri - 1^; (10) 

This allows to define the life time r, 77- = The notation << ... >> stands for the white 

noise average and A is the momentum cut-off. The real part of the self energy diverges in 
the hmit u — 0. The derivative of the self energy introduces the wave function renormal- 
ization Z"^ (or the quasi peirticle weight): Z~^{uj,X;A) = [1 — dcjT,ii{u), fj,, X, A; p)\cj=o\ 
As a result we replace the chemical potential by the renormalized chemical potentiah/Lt — > 
/X = Z[iJ, — Sr(0, fi, A, A;p))] Similarly the velocity v is replaced by : f — >■ vZ{uj, A; A) The 
averaged Green's functions for particles GP{k; cu) and anti-particles D"-{k; u) is given by : 
G^{k; u) = ds&{k; s = t-t')e'^'+J^^, ds&{k; s = t-t')e"^' = &{k; u, +)+&{k; u, -) 
and W{k; cu) = dsG^{k; s ^t- t')e'^' + dsG^{k] s ^t- t')e'^' = G^{k; oj, +) + 
G°'{k; cu, —). The averaged Green's functions can be expressed in terms of the spectral func- 
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tions A± for particles and B± for antiparticles . 

G^^k- ^, +) ^ r dz- -; G^ik; u, -) ^ r dz- 



a) — — fi) + ie J-oo Cj — {E{k) — ji) — ie 

A^{k,co) = (^) -J'-''-^^^^ ;A4k,u) ' ^ 

+ ^ ' ^ MtttVa ^z?/'^,^ ..y\2 i M ^2' ^ ' 
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oo 



G^(A;;c^,+)= / dz- -.G''^k;U,-)= r dz- ^'^^'^^ 



oo UJ + {E{k) + jj,) + ie J~oo Cj + {E{k) + jj,) - ie 

B+{k,u) = (7—) ,. , ,^,r^ , ^^2 , /M2 '^-(^'^ " 



AnT\cJ + {E{k) + ^^)f + {j-y Ht^t^Cj + {E{k) + + {If' 

(11) 



V-Computation of the current in the averaged single particle approximation 

We replace the ground state |0 > with the average effective ground state |0o > charac- 
terized by the spectral functions A± and B±. We construct the evolution operator due to 
the external potential H'^^^{t) which acts on the effective ground state |0o >• We use the 
interaction picture and compute the induced current 5Ji{r,t) to linear order in the vector 
potential A{t) (see eq.l): The Hamiltonian which describe the coupling of light to matter is 
given by : 

H'-\t) = I d'r[J,if,t)A,{f,t) + Mf,t)A,if,t)] (12) 

We use the interaction picture and compute the induced current SJi{f,t) by the vector 
potential A{t) (see eq.l): 

SJi{f,t) = ^ f dt' < 0o|[Ji(r,t),//^^*(t')]l0o > (13) 

The current operator is defined by the variation of the Hamiltonian with respect the vector 
potential: 

oAi{r,t) 



oA2{r,t) 



(14) 



we obtain the linear response for the conductivity ai^i{q,Q) |l4|. 
cxM(g,^^) = ^^ 
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/oo 
-oo 



= t - = Y^[t - t'] < <^o|[ii(g,t), Ji(-g,0]Po > 



(15) 



^[t — t'] is the step function which is one for t > t'. The current operator Ji(r) is build 
from the four components Jf'^(g; t), Jf'"(g; i), Jf'^(g; f) Jf '"(g; i) and the spinor matrix el- 
ements : Wf'^{k + q,k) =< u^+\k + q)\a^\u^+\k) >;Wt\-k - q, -k) =< u^-\-k - 
q)\a^\u^-\-k) > ■,W^'''{-k - q,k) =< u^-^{-k - q)\a^\u^+\-k) >; Wr{k + q,-k) =< 
M(+)(fc + g)|(7>(-)(-fc) > 

jrm = J ^,c\k+q,t)wr{k+q,k)c{it) 

•^"'"^^"'^^^ / -^2B{-k-q,t)Wr{-K-l-k)B\-k,t) 
Jrm^ J -0^^C^{k + q,t)Wr(k + q-k)B^(-k,t) 
Jrm^ I ^^B{-k-q,t)Wr{-k-q,k)C{k,t) 



(16) 



We can express the conductivity in terms of the four currents build from the particles p and 
anti-particles a :Jf'''(q;t) ,J?'''{q;t) Jr(q;t), J^'''{q;t). 

Ri^i {q, n;p,p; p, p) + Ri^i {q, n;p,a; a, p) + Ri^i {q,n;a,p; p, a) + Ri^i {q,n;a,a;a,a) 



in 

(17) 



Using the explicit form of the spectral functions A±,B± given in equation (10) we obtain 
the conductivity. 

a) The conductivity in the limit of infinitesimal disorder with a vanishing 
chemical potential 

The only finite contributions are given by the combination B_A^ , the other two spectral 
functions B. and A_ are zero. 



<7i,i(!J; l'^0) = —J_^tLjJ^ EdEl 
2h J-oo Jo il 



B_{E,Lj)A+{E,LJ + n) - B_{E,LJ + n)A+{E,u)]. 
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-^[/F.i?.(-^)(1 - /f.d.(^)) + /f.d.(^)(1 - /f.d.(-^)]|t-^o ^ 

(18) 

b) -The conductivity for a finite chemical potential > in the limit of in- 
finitesimal disorder 

For this case = and we have three nonzero spectral functions ,A_^.,A_ and B_. 

2h J-oo Jo il 
-B_ {E, u + Q)A+{E, u) + B_ {E, uj)A+{E, uj + fl) 

(19) 

2 

We find that the the conductance at T = is given by the universal value ^ for Q > 2jj,. 
For a finite chemical potential the metallic behavior is given by ^\fi\S{il). In figure 1 we 
show the conductivity for the entire range of frequencies for the case that the elastic mean 
free path is large. 

c) - The conductivity in the limit of finite disorder with the finite life time 

case lei. — VT at a low frequency Q < 

In the limit of in the limit — )■ we can limit ourself only to the conductivity of the 
conduction band which is given by given by i?i,i(0, fl;p,p;p,p). We find Drude like behavior 
given by figure 2. 

Where kp is given by kp = We observe that the conductivity contains the factor half. 
The origin of the factor half is due to the angle dependent vertex cos^[x(A;)]. In addition we 
remark that Ladder correction will replace the scattering time r and therefore the elastic 
mean free path with the transport time Tfr > r and therefore I with If^ > I- In figure 2 we 
have plotted the conductivity for this case. 
Conclusion 
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FIG. 1: The conductivity for the entire frequency 

To conclude we have computed the optical conductivity for the Weyl Hamiltonian which 
describe the surface excitations of the TI for the entire range of frequencies. Due to the time 
reversal invariance d"^ = —1 the backscattering is prohibited and the model belongs to the 
sympectic ensemble justifying the use of the single particle approximation. The universal 

2 

value of the conductivity is given by a{u > fi) = For a finite chemical potential and 
finite disorder we find that the conductivity is given hj a{Q < 2 



^'\) = kwpTT- In the limit 
of ~ self consistent calculations performed by one of us 17| show that the conductivity 
is given by a{n ~ 0, f2 = 0) 
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